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additional  effects  of  flow  velocity,  ripple  length,  and  deviations  of  the 
ripple  profile  from  normal  equilibrium.  Average  rates  of  energy  dissipation 
are  calculated.  Results  show  f  to  be  a  complicated  function  of  0  and 
other  parameters.  In  particular,  the  instantaneous  stress  is  not  simply 
related  to  the  instantaneous  velocity.  Some  salient  features  of  f(0)  are 
described,  qualitatively,  by  a  simple  model. 
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MEASUREMENTS  OF  OSCILLATORY  DRAG  ON  SAND  RIPPLES 


Karl  E.  B.  Lofquist 
National  Bureau  of  Standards 
Washington ,  D.C. 


ABSTRACT 

Measurements  have  been  made  of  drag  on  naturally  rippled  sand  beds 
in  an  oscillatory-flow  water  tunnel.  A  partition  splits  the  tunnel  into 
two  parallel  channels  with  equal  cross  sections  and  volume  rates  of 
flow,  one  with  a  smooth  flat  rigid  bottom  and  the  other  contai  Ing  the 
sand  bed.  Roughly,  the  difference  in  the  two  bottom  drags,  the  one 
known,  is  equal  to  the  difference  in  the  net  pressure  forces  on  the  water 
in  the  two  channels,  which  is  obtained  from  measurements  of  pressure 
differences  across  the  partition  at  each  end  of  the  sand  bed.  Each 
experiment  provides  the  drag,  or  average  bottom  stress,  on  the  rippled 
sand  bed  as  a  function  of  the  phase,  6,  of  the  sinusoidal  flow.  From  a 
first  set  of  thirteen  experiments  with  a  medium  sand,  a  stress  coeffi¬ 
cient,  f(0),  is  presented  in  three  families  of  curves  which  explore  the 
additional  effects  of  flow  velocity,  ripple  length,  and  deviations  of 
the  ripple  profile  from  normal  equilibrium.  Average  rates  of  energy 
dissipation  are  calculated.  Results  show  f  to  be  a  complicated  function 
of  0  and  other  parameters.  In  particular,  the  instantaneous  stress  is 
not  simply  related  to  the  instantaneous  velocity.  Some  salient  features 
of  f (0)  are  described,  qualitatively,  by  a  simple  model. 


I.  INTRODUCTION 

1.  General 

Water  waves  on  the  sea  surface  and  sand  ripples  on  the  sea  bed 
interact  by  processes  of  drag.  Surface  waves  provide  an  oscillatory 
flow  above  the  bed  which  produces  and  maintains  bed  forms  by  local 
movements  of  sand.  The  sand  is  scoured  from  the  bed  by  tangential 
stresses  under  flows  which  include  separation  and  vortex  formation  and 
which  largely  are  shaped  by  the  pressure  field  around  the  bed  form. 

This  field  produces  a  distribution  of  normal  stress  over  the  surface 
of  the  bed.  The  horizontal  force  of  these  tangential  and  normal 
stresses  on  a  large  area  of  the  bed,  divided  by  that  area,  defines  the 
average  bottom  stress.  This  is  a  function  of  time  which  becomes  perio¬ 
dic  as  the  bed  forms  become  stable.  When  the  bed  form  is  two  dimen¬ 
sional  and  periodic  in  space,  the  bottom  stress  is  obtained  by  an 
average  over  a  single  ripple  length.  The  reverse  effect  of  the  ripples 
on  the  surface  waves  is  less  direct  and  relatively  minor,  but  is  signi¬ 
ficant.  The  bottom  stress  combines  with  the  oscillatory  flow  above  the 
bed  to  dissipate  energy.  Ultimately,  this  energy  is  supplied  by  the 
surface  waves  which  decay  in  the  process.  Also,  the  bottom  stress  Is  a 
factor  in  the  refraction  of  surface  waves.  This  paper  describes  labora¬ 
tory  measurements  of  bottom  stress  on  naturally  rippled  sand  beds  as 
function  of  time. 
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2.  Definition  of  Terms 

Terms  introduced  in  the  preceding  paragraph  with  others  to  be  used 
in  this  study  are  defined  by  a  few  equations.  Above  the  bed  and  beyond 
its  effects  the  velocity  of  flow  is 

u»  =  (^P)  sin  (^P)  *  u  sin  9  (X) 

where  t  is  the  time,  T  is  the  period  and  2a  is  the  orbital  diameter  of 
the  water  motion.  Equation  (1)  defines  the  maximum  velocity,  U,  and  the 
phase,  6.  The  bed  surface  profile,  assumed  two  dimensional  and  periodic 
with  length  X,  has  elevation  h(x,0),  where  x  is  the  horizontal  coor¬ 
dinate.  The  average  bottom  stress  is  then  defined  by 

X 

T(e>  -  x  |  f  th  +  (Ph-P..>!Jj  I  d*  (2) 

0 

where  is  the  local  tangential  stress  on  the  bed,  and  (p^-p^)  is  the 
difference  between  pressures  at  the  bed  surface  and  at  any  fixed  eleva¬ 
tion  beyond  the  influence  of  the  bed._  Separate  integrals  for  the  two 
terms  in  the  integrand  would  resolve  t(9)  into  components  expressing  the 
effects  of  tangential  and  normal  stresses,  the  latter  being  a  form  drag. 
However,  such  a  separation  is  generally  impossible.  A  stress  coeffi¬ 
cient  , 

<» 

where  p  is  the  density  of  the  water,  has  been  made  dimensionless  by 
dividing  by  pU2 ,  rather  than  by  pu2(6),  in  order  to  keep  f(8)  propor¬ 
tional  to  t(0)  throughout  the  cycle.  The  maximum  of  f(0)  is  denoted 
f  .  The  energy  dissipated  per  unit  area  of  bed  per  period,  or  twice 
t¥iat  per  half  period,  is  given  by 

T/2 

AE  *  2  f  t  u#  dt  (4) 

J  0 

Then,  inserting  equations  (1)  and  (3)  into  equation  (4)  and  defining 
the  average  coefficient, 

fl  *  ~  f  f(8)  sin  0  d0  (5) 

J  o 

produces. 


Thus,  the  rate  of  energy  dissipation  is  proportional  to  f  j . 


(6) 
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3.  Past  Work 


In  the  past,  f(0),  f  and  f1  have  been  measured  or  estimated  by 
various  methods.  These  methods  are  here  divided  into  two  groups  accord¬ 
ing  to  whether  the  surface  of  the  bed  was  (1)  artificially  rigid  or 
(2)  naturally  rippled. 

(1)  With  rigid  surfaces:  Bagnold  (1946)  calculated  fj  by  the  work 
needed  to  oscillate  a  metal  cusp-shaped  profile  in  still  water. 

Kalkanis  (1964)  oscillated  roughened  plates  in  still  water  and  observed 
flows  which  could  be  analyzed  to  provide  f(0).  Using  oscillatory-flow 
water  tunnels,  Jonsson  and  Carlsen  (1976),  in  two  experiments,  measured 
velocity  profiles  over  a  flat  roughened  bed  to  obtain  f(0),  while 
Kamphuis  (1975)  obtained  f  on  a  roughened  shear  plate  in  terms  of 
relative  grain  roughness  and  a  Reynolds  number. 

(2)  With  naturally  rippled  beds,  fj  only  has  been  obtained.  By 
observation  of  the  decay  of  surface  waves,  fj  has  been  calculated,  in 
wave  tanks,  by  Savage  (1953)  and  by  Inman  and  Bowen  (1962)  and,  in  field 
studies,  by  Bretschneider  (1954)  and  by  Iwagaki  and  Kakinuma  (1967). 
Values  of  fj  from  these  field  studies  are  summarized  by  Vitale  (1979) 
and  vary  widely.  Cars tens,  Neilson  and  Altinbilek  (1969)  measured  the 
energy  input  required  to  maintain  flows  over  sand  beds  in  an  oscillatory- 
flow  tunnel  restricted  to  a_single  period  of  3.56  seconds  and,  within 
this  limitation,  presented  f1  as  function  of  grain  size,  D,  and  of  a/D. 
Vitale  (1979)  has  further  analyzed  and  interpreted  their  results. 

In  sum,  except  with  rigid  flat  roughened  beds,  measurements  of 
bottom  drag  remain  fragmentary.  In  particular,  no  previous  measurements 
of  f ( 0)  with  naturally  rippled  beds  have  been  reported. 

4.  Scope  of  the  Present  Study 

To  help  fill  this  gap,  the  present  study  measures  f(0)  in  labora¬ 
tory  experiments  with  naturally  rippled  sand  beds.  Experiments  are 
still  in  progress,  but  results  from  thirteen  o^f  them  are  presented,  and 
preliminary  conclusions  are  drawn.  Values  of  fj  are  calculated,  and  some 
features  of  f(0)  are  discussed  in  terms  of  a  simple  model. 


II.  EQUIPMENT  AND  PROCEDURES 
1.  The  Tunnel  and  its  Modifications 


The  experiments  are  carried  out  in  an  oscillatory-flow  water 
tunnel  described  by  Lofquist  (1977).  Briefly,  it  is  of  U-tube  design, 
and  two  tight  fitting  pistons  in  cylinders,  side  by  side  at  one  end, 
moving  sinusoidally  in  unison,  drive  water  back  and  forth  through  a 
horizontal  test  section  and  into  and  out  of  two  reservoirs  open  to  the 
air  at  the  other  end.  The  test  section  is  255  cm  long,  21  cm  wide,  and 
has  a  mean  depth  of  flow  of  30  cm.  The  orbital  diameter,  2a,  and  the 
period,  T,  can  be  selected  arbitrarily  below  250  cm  and  above,  approxi¬ 
mately,  1.5  seconds.  Ranges  in  the  present  experiments  are: 

34.8  cm  <  2a  <  76.6  cm,  and  2.70  s  <  T  <  5.93  s. 


I 


COASTAL  ENGINEERING- 1980 


reservoirs  pistons 


tops  partition  taps 


gap  m  partition 


Figure  1.  Schematic  drawings  of  the  tunnel  and  test  section, 

(a)  side  view,  (b)  top  view,  (c)  details  of  the  two  ends 
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The  tunnel  has  been  modified  for  the  present  study  as  shown  in  the 
three  schematic  drawings  of  Figure  1(a) , (b) ,  (c) .  Of  these,  only  (c)  is 
drawn  to  scale.  The  channel  has  been  split  into  two  channels,  each 
10.1  cm  vide,  by  a  vertical  partition  (Fig.l,b).  This  partition  sepa¬ 
rates  each  side  from  the  other  all  the  way  from  the  pistons,  at  the 
right,  to  a  point  beyond  the  left  end  of  the  test  section  where  the 
partition  ends  at  a  screen  and  there  is  a  bypass  over  the  channel  top 
(Fig. l,c) .  (A  screened  bypass  enlargement  was  cut  in  the  partition  just 
before  the  last  three  experiments.)  One  channel  has  a  flat  smooth  metal 
bottom,  while  the  other  contains  the  rippled  sand  bed  (Fig.l.b).  In 
both  channels,  beyond  the  test  section  at  either  end,  pressure  taps  have 
been  mounted  in  the  channel  walls.  These  are  indicated  in  all  the  draw¬ 
ings.  A  pitot  tube  has  also  been  mounted  which  can  be  raised  and  lowered 
and  which,  with  oppositely  directed  openings  at  its  tip  serving  alter¬ 
nately  as  static  and  dynamic  holes,  can  measure  continuously  the  velocity 
of  the  reversing  flow.  The  pressure  taps  in  the  walls  and  the  pitot  tube 
openings  are  connected  to  two  differential  pressure  transducers  which  are 
located  as  indicated  in  Figure  l,(a).  The  connecting  conduits  are  not 
shown.  The  purpose  of  these  modifications  will  appear  in  the  derivation 
of  the  stress  measurement  in  section  11,3  below. 

2 .  End  Conditions,  Natural  and  Strained  Profiles 

The  analysis  in  section  11,3  requires  that  the  areas  of  the  four 
cross  sections  at  the  pressure  taps  be  equal,  so  that  migrations  of  the 
sand  profiles  into  the  conduits  must  be  prevented.  Therefore,  rigid 
crest-shaped  objects  have  been  mounted  at  each  end  of  the  sand  bed  and, 
also,  just  before  the  first  screens  in  order  to  stabilize  the  profiles 
and  to  keep  the  floor  below  the  taps  free  of  sand.  These  objects  are 
shown  in  Figure  l,(c).  The  rigid  crests  at  the  ends  of  the  bed  are 
either  of  two  sizes,  selected  best  to  match  the  sand  crests  between 
them.  These  rigid  end  crests  condition  the  profiles  since  the  number  of 
sand  ripple  crests  between  them  is  some  integer,  M,  and  the  fixed  dis¬ 
tance  between  them  must  equal  A(M+1).  Studies  have  shown  (Lofquist, 

1978,  Fig.  22  and,  more  fully,  Brebner,  elsewhere  in  these  proceedings) 

that  an  unconstrained  two-dimensional  profile  adjusts  A /a  to  a  natural 
value  (X/a)  which  is  a  function  of  the  degree  of  sand  agitation  and 

grain  size  But  for  low  and  moderate  agitation  is  approximately  equal  to 

1.33.  Thus,  in  the  tunnel,  with  X  and  a  selected  such  that  X/a  -  (A/a)^ 
the  rigid  end  crests  create  minimum  distortion  of  the  sand  profile  between 
them. 


With  X  fixed  by  a  selected  number  of  sand  ripple  crests  between  the 
fixed  end  crests,  it  is  possible,  by  varying  a,  deliberately  to  force 
X/a  away  from  (X/a)  and  so  to  produce  a  "strained"  profile.  Such  pro¬ 
files  are  described1^  Lofquist  (1978),  V,5.  Briefly,  as  X/a  is  reduced 
below  1.33,  crests  steepen  and  the  profile  appears  "compressed"  until  at 
around  X/a  *  1.1  somewhere  a  crest  shrinks  and  disappears,  and  X  abruptly 
increases.  As  X/a  is  Increased  above  1.33,  crests  flatten  and  the  pro¬ 
file  appears  "extended",  with  low  humps  forming  in  the  troughs,  until 
at  around  X/a  -  1.6  one  of  these  hraps  rises  to  become  a  new  crest,  and 
X  abruptly  decreases.  The  deliberate  variation  of  X/a  permits  an 
investigation  of  the  effect  of  profile  strain  on  f(6)  and  whether  a 
natural  profile  gives  f(0)  any  special  character. 
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3.  Derivation  of  the  Measurement  of  Stress 


A  simple  analysis  provides  the  average  stress  on  the  sand  bed  in 
terms  of  the  pressures  at  the  four  taps  and  the  velocity  profiles  mea¬ 
sured  by  the  pitot  tube.  Symbols  used  in  the  analysis  are  identified  as 
applying  either  to  the  channel  with  the  sand  bed  or  with  the  smooth 
bottom  by  a  subscript  "s,f  or  "o"  respectively.  These  symbols  are:  L, 

H  and  b,  the  length,  depth  and  width  of  the  channels  (the_same  for  each); 
p(0)  and  p(L) ,  the  pressures  at  the  taps  at  x  ■  0  and  L;  t  and  t  ,  the 
average  stresses  on  the  bottom  and  on_the  smooth  perimeter  of  thl  walls 
and  top;  Q,  the  volume  rate  of  flow;  u2 ,  the  average  of  u2  over  a  cross 
section,  where  u  is  the  local  velocity. 

Then,  by  momentum  balance  in  each  of  the  channels, 

bH  [p(0)-p(L>]  -  bL  7  -  L(2H  +  b)  7  -  pLQ  +  pbH  [72(L)-72(0)]  (7) 

s  s  ps  s  s 

and 

bH  [p(0)-p(L))o  -  bL  7o  -  L(2H+  b)  tpo  *  pLQo  +  pbH  [72(L)-u2(0)  ]q  (8) 


where  a  dot  denotes  a  derivative  with  respect  to  time.  Since  Q  is  known 
and  t  can  be  estimated,  t  can,  in  principle,  be  determined  by 
equation  (7).  This  approach  was  attempted  before  the  channel  was  divided. 
But  each  side  of  equation  (7)  is  dominated  by  its  first  term,  and  Lx 
is  easily  obscured  by  noise  and  scatter  in  measurements  of 
H[p(0)-p(L)]  .  #  Instead,  equation  (7)  is  subtracted  from  equation  (8). 

In  doing  so,  Q  -  Q  vanishes,  since  Q  and  Q  are  made  equal  by  con¬ 
struction  of  the  tunnel  with  two  pistons  of  equal  area.  Also,  t  -  t 
can  be  discarded,  since  the  terms  are  of  nearly  equal  magnitude  End  eaEfi, 
as  estimated  by  laminar  theory,  is  small  compared  with  x  as  eventually 
determined.  Finally,  since  the  development  of  the  velocity  profile  over 
the  smooth  bottom  is  relatively  small,  the  last  term  in  equation  (8)  can 
be  discarded.  The  subtraction  then  provides. 


Ts  ■ 


L  [(VPs>L  "  (po_ps)01 


*£  [7|(l)-?(0)] 


(9) 


where  the  pressure  differences  are  between  pairs  of  taps  at  x  *  L,0. 

With  x  calculated  by  laminar  theory,  results  of  the  experiments 
eventually  show  that  x  dominates  the  left  side  of  equation  (9)  and  that 
the  right  side  is  dominated  by  the  first  pressure  difference,  with  the 
second  kept  relatively  small  by  the  action  of  the  by-pass.  The  last 
term,  which,  like  x  ,  proves  to  be  small  compared  with  x  ,  is  provided 
by  the  single  pitot°tube  at  x *  L  with  the  assumption  thaf,  with  z  the 
vertical  coordinate,  us(x,z,0)  satisfies  the  condition  of  symmetry, 


u  (Ofz,0)  *  -  u  (L,z,6-Hr) 

S  8 


(10) 
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4.  Pressure  Measurement 


A  differential  pressure  transducer  mounted  over  the  test  section 
near  the  taps  at  x  ■  L  (Fig.l,a)  is  alternately  connected  with  the  two 
pairs  of  taps  at  x  *  0,L  by  suitable  conduits  and  valves  (not  shown  in 
the  drawings).  Next  to  it  is  mounted  a  second  pressure  transducer  con¬ 
nected  with  the  pitot  tube.  Signals  from  the  transducers  are  recorded 
on  paper  tape  as  functions  of  time.  Pen  deflections,  in  mm,  are  trans¬ 
lated  into  pressure  units  by  a  calibration  in  which,  immediately  follow¬ 
ing  an  experiment,  a  simple  device  provides  selected  hydrostatic  pressure 
differences  to  a  transducer  while  the  recorder  runs  with  the  same  ampli¬ 
fication  as  before. 

Filters  are  required  to  reduce  noise  in  the  transducer  output 
signals.  These  filters  consist  of  six  sand-filled  glass  loops  one  of 
which  is  included  in  each  of  the  conduits  between  the  transducers  and 
the  four  pressure  taps  and  the  two  pitot  tube  openings.  This  filtering 
reduces  the  amplitude  and  slows  the  response  of  the  output  signal. 

These  distortions  must  be  removed  by  calibration  of  another  type.  With 
the  pressure  at  one  tap,  or  pitot  opening,  fixed  by  draining  the  water 
in  the  tank  to  the  level  of  an  open  hatch,  the  other  is  connected  by  a 
flexible  tube  to  a  vessel  which  is  filled  with  water  to  this  same  level. 
The  vessel  is  then  oscillated  vertically  in  simple  harmonic  motion  with 
selected  amplitude,  h  ,  and  period,  T  ,  to  supply  a  known  (non  hydro¬ 
static)  sinusoidal  pressure  difference  to  the  transducer.  Denoting,  for 
given  h  ,  the  ratio  of  the  amplitude  of  the  output  signal  as  T  -►  ®  to  that 
at  by  and  the  phase  lag  by  <f>,  calibrations  are  obtained  Sn  the  form, 

T  2  ($2-1)  *  F,(T  )  (11) 

c  1  c 

and 

Tc  tan  *  -  F2(Tc)  (12) 


where  F1  and  F2  are  independent  of  h  and  are  slowly  varying  functions  of 
T  or  constant.  Although  subject  to  adjustments  and  slight  drift  over 
long  periods  of  time,  representative  constant  average  values  for  both 
pairs  of  taps  are,  Fj  *  A. 3  sec2,  F2  -  2.5  sec,  and  for  the  pitot  tube, 

Fj  a  9.7  sec2,  F2  =  4.3  sec. 

5.  Processing  the  Data 

As  the  first  term  on  the  right  side  of  equation  (9)  has  been  found 
to  dominate  the  last  term,  the  processing  of  records  from  the  pressure 
taps  has  differed  from  that  of  records  from  the  pitot  tube  and  will  be 
considered  first.  The  primary  records  from  the  pairs  of  taps  at  x  *  0,L 
are  the  pen  deflections,  in  mm,  which  are  denoted  P(0,6)  and  P(L,6)  respec¬ 
tively.  Figure  2  shows  P(L,0)  and  P(O,0)  for  a  typical  experiment  (no.  2 
in  Table  1  below).  In  each  case  residual  noise  has  been  averaged  out  to 
produce  a  smoother  curve.  The  deflections  of  P(O,0)  are  kept  relatively 
small  by  the  action  of  the  by-pass.  The  phase  is  marked  for  0  -  0,  tt, 

2 it,  A  digitizing  pad,  or  table,  reads  the  analog  tape  records  at  each 
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0  ▼  2  ir 


6 


Figure  2.  A  typical  record  of  transducer  output 
from  the  two  pairs  of  pressure  taps. 


10  degrees  of  phase  and  stores  the  data  In  a  microcomputer.  The  computer 
then  processes  the  data  In  the  following  steps: 

(1)  As  required  in  equation  (9),  the  difference  in  the  records  is 
taken  to  fora 

0(6)  -  P  (L ,  0 )  -  P(O,0)  (13) 


(2)  Ideally,  the  stress  at  any  time  should  be  the  negative  of  that 
a  half  cycle  earlier  or  later  in  the  opposite  direction »  so  that  0(6) 
should  equal  -  D(0  +  *).  This  condition  is  only  approximately  realized. 
Therefore,  D(0)  is  averaged  with  -  D(0  ♦  it)  to  produce  the  function 

ff(e)  -  \  [D(e)  -  D(e  +  «)i  (U) 

which  changes  only  in  sign  when  advanced  a  half  cycle.  In  Figure  3 
below,  for  each  of  the  thirteen  experiments,  points  giving  D(0)  and 
-  D(0  +  k)  are  plotted  at  each  10  degrees  of  phase,  and  D(6)  is  shown. 
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(3)  D(6)  is  resolved  into  a  sum  of  sinusoidal  components  by  Fourier 
analysis.  Both  sine  and  cosine  terms  are  present  with  arguments  n0 , 
where  n,  an  integer,  is  the  number  of  the  harmonic,  but  since  D(0)  - 

-  D(8  +  it),  even  harmonics  are  absent,  and  n  is  odd. 

(4)  Successively  taking  T  *  T/n,  the  calibration  corrections  for 
the  filter  distortions  given  by  equations  (11), (12)  are  applied  to  each 
component ,  and 

(5)  The  series  is  reassembled  to  provide  a  corrected  D(0).  This 
reassembly  has  included  terms  as  far  as  the  7th  harmonic  while  checking 
that  the  increment  from  the  9th  harmonic  is  insignificant. 

(6)  By  use  of  the  hydrostatic  pressure  calibrations,  D(0)  is 
expressed  in  pressure  units  to  provide  the  first  term  on  the  right 
side  of  equation  (9) . 

The  records  from  the  pitot  tube,  also  contained  on  the  paper  tapes, 
have  been  processed  manually,  without  computer  or  Fourier  analysis.  The 
last  term  in  equation  (9)  has  been  approximated  using  overall  filter 
corrections,  those  provided  by  equations  (11)  and  (12)  for  n  *  1.  Errors 
associated  with  this  approximation  cannot  be  serious,  since  this  term  has 
been  found  to  be  small,  its  range  less  than  8%,  and  typically  less  than 
5%,  of  the  main  term. 

The  remaining  unknown  term  in  equation  (9) ,  t  ,  is  calculated  by 
the  theory  for  laminar  oscillatory  flow  over  a  plane  surface.  Its  range 

has  proven  to  be  less  than  3%  of  that  of  t  . 

s 


IV.  RESULTS 


1.  Summary  of  Experiments 

The  major  results  of  this  paper  are  curves  of  f  as  function  of  0  and 
of  the  parameters  U,  A,  A /a  from  the  thirteen  experiments  analyzed  so  far. 
In  Table  1  these  experiments  are  Identified  by  number,  in  chronological 
sequence,  and  by  values  of  their  parameters  M,  A,  a,  U,  and  including  the 
redundant  T  -  2xa/U.  The  trough  to  crest  ripple  height,  n,  has  also  been 
added  (except  where  the  record  remains  on  undeveloped  film) .  H  is  the 
number  of  sand  crests  between  the  fixed  end  crests  and  determines  A,  as 
described  in  section  11,2.  Experiments  for  which  A/a  differs  from 
(A/a)  ,  that  is,  those  with  strained  profiles  (sec. II, 2),  are  denoted  by 
asterisks.  Experiment  "C"  is  a  composite  of  nos.  2  and  26  as  described 
in  section  IV, 2.  All  thirteen  experiments  have  used  quartz  sand  with 
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Table  1.  Parameters  for  the  Thirteen  Experiments 
*  0.55  mm 


Exp.  no. 

M 

A,  cm 

a,  cm 

U,em/s 

T,s 

n,cm 

1 

7 

31.8 

23.9 

27.0 

5.57 

7.3 

2 

7 

31.8 

23.9 

40.6 

3.71 

6.9 

3 

7 

31.8 

23.9 

33.1 

4.55 

7.8 

4* 

7 

31.8 

25.5 

33.3 

4.83 

7.2 

5* 

7 

31.8 

21.9 

33.1 

4.16 

6.7 

6* 

7 

31.8 

20.6 

33.15 

3.90 

5.7 

7* 

7 

31.8 

27.7 

33.2 

5.25 

7.6 

13 

10 

23.1 

17.4 

40.5 

2.70 

4.8 

18 

5 

43.5 

31.9 

40.6 

4.94 

9.3 

24 

4 

52.5 

38.3 

40.6 

5.93 

11.9 

26 

7 

31.8 

23.9 

40.5 

3.72 

- 

27 

7 

31.8 

25.9 

52.1 

3.13 

- 

28 

7 

31.8 

28.9 

63.6 

2.86 

- 

C 

7 

31.8 

23.9 

40.5 

3.71 

- 

*With  strained  profiles,  A/a  j  (A/a)n 


2.  Reliability 

Interpretation  of  the  results  to  follow  requires  some  assessment  of 
their  reliability.  The  most  convenient  measure  of  reliability  is  the 
extent  to  which  the  stress,  and  related  measurements,  change  only  in 
sign  as  Q  is  increased  by  it,  that  is,  the  agreement  between  D(e)  and 
-  D(0  +  it)  defined  in  sec.  11,5  and  by  equation  (13).  In  Figure  3, 

D(9) ,  -  D(0  +  it)  and  their  average  D(0)  (eq.  14)  are  shown  for  each 
experiment,  identified  by  its  number  in  Table  1.  The  scales  of  D(e) 
and  -  D ( 0  +  tt)  ,  are  in  tape  units  before  calibration  and  are  of  no 
interest  beyond  the  indicated  zero  marks.  Discrepancies  between  D(o) 
and  -  D( 0  +  tt)  remain  within  10%  for  the  most  part  but  are  occasionally 
considerably  larger.  Their  cause  is  unknown.  They  are  most  plausibly 
associated  with  pressure  buildups  behind  the  screens  and  resulting 
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flows  through  the  by-pass.  Pressure  buildups  are  caused  especially  by 
sand  in  suspension  being  carried  to  and  accumulating  against  a  screen 
during  a  part  of  the  cycle  and  become  large  enough  to  require  limitations 
on  U  in  order  to  forestall  damage  to  the  screens.  The  discrepancies  may, 
in  part,  also  derive  from  slight  asymmetries  in  the  sand  profiles  from 
one  end  of  the  bed  to  the  other,  particularly  near  the  fixed  end  crests. 
Whatever  their  cause,  such  discrepancies  reduce  confidence  in  an  experi¬ 
ment  and  if  too  large  may  cause  it  to  be  discarded. 

Three  of  the  experiments,  nos.  26,27,28,  were  performed  (after  the 
oral  presentation)  following  an  enlargement  of  the  by-pass  to  reduce  the 
pressure  difference  (P0~Ps)q  9)  especially  at  higher  values  of  U. 

Except  for  this  enlargement,  exps.  2  and  26  are  ostensibly  the  same,  and 
have  been  combined  into  a  composite  experiment  numbered  "CM.  Figure  4 
shows  f(0)  for  exps.  2,  26,  and  "C" .  These  curves  agree  reasonably  well 
over  most  of  the  range  but  diverge  substantially  near  times  of  flow 
reversal.  Exps.  27  and  28  replace  two  experiments,  included  in  the  oral 
presentation,  which  had  the  same  values  of  U,  but  had  A/a  *  1.33,  rather 
than  (A /a)  as  in  exps.  27  and  28.  At  these  high  values  of  U,  (A/a)^  is 
significantly  less  than  1.33.  In  these  two  earlier  experiments,  ranges 
of  P(0,6)  exceeded  one  half  those  of  P(L,0)  and  discrepancies  between 
D(0)  and  -  D(0  +  tt)  were  substantially  larger  than  in  exps.  27  and  28, 
and,  perhaps  in  consequence,  their  curves  of  f ( 0 )  were  quite  different. 


0 


w 


8 


Figure  4.  f(0)  for  experiments  2  and  26  and  their  composite  "C" . 
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3.  The  Curves  of  f(e) 

The  curves  of  f (0 )  from  the  twelve  experiments  (after  combining 
nos.  2  and  26),  for  selected  parameters,  are  shown  In  Figures  5,6,7.  In 
Figures  5  and  6  the  profiles  are  natural,  with  A/a  *  (A/a)  ,  and  in 
Figure  7  they  are  strained.  (A/a)  ■  1.33  in  all  experiments  except 
nos.  27  and  28,  in  Figure  5,  where  it  is  reduced  by  the  large  values  of 
U.  In  each  figure,  or  family,  the  curves  are  identified  at  their  left 
by  their  numbers  in  Table  1  and  are  arranged  in  order  of  monotonic 
increase  or  decrease  in  the  defining  parameter.  In  Figures  5  and  6  this 
parameter  is  U  and  A,  respectively,  and  in  Figure  7  it  is  A/a.  The 
numerical  value  of  the  varying  parameter  is  given  at  the  right  of  each 
curve,  and  conditions  and  values  of  parameters  held  constant  are  shown 
above  each  family  or  figure.  Although  f(0)  is  determined  by  its  behavior 
over  the  restricted  range  o  <  0  <  it  ,  since  f(0  +  tt)  -  -  f(0),  it  is  here 
presented  over  the  entire  cycle  -  ir/2  <  0  <  3x/2  so  that  features  of  f(0), 
such  as  maxima  and  minima,  can  more  easily  be  followed  from  one  curve  to 
the  next  as  the  parameter  is  varied. 

4.  Values  of  f,  and  f 

- 1 - m 

With  equation  (5),  values  of  the  measure  of  energy  dissipation, 
have  been  calculated  from  the  twelve  curves  of  f(0)  already  presented. 
Table  2  lists  these  values  in  the  same  order  as  the  curves  in  Figures  5, 

6,7.  For  comparison,  maximum  values  of  f(0),  f  ,  have  been  included. 

m 


Table  2.  Values  of  f^,  proportional  to  the  rate  of  energy 

dissipation,  and  f  for  the  curves  in  Figures  5,6,7. 


Natural  profiles 

Strained  profiles 

A 

-  31.8 

cm 

U  « 

40.6 

cm/s 

A 

*  31.8 

cm 

U  ~ 

33.2 

cm/s 

U  ,cm/s 

f. 

f 

A  ,cm 

f. 

f 

X/a 

f. 

f 

1 

m 

1 

m 

1 

m 

27.0 

.122 

.299 

52.2 

.129 

.282 

1.15 

.113 

.253 

33.1 

.102 

.236 

43.5 

.115 

.259 

1.25 

.105 

.228 

40.5 

.083 

.180 

31.8 

.083 

.180 

1.33 

.102 

.236 

52.1 

.075 

.218 

23.1 

.115 

.452 

1.45 

.107 

.228 

63.6 

.072 

.348 

1.54 

.086 

.227 

The  combination  U  -  40.6  cm/s  and  a  *  23.0  cm  occurs  In  the  observa¬ 
tions  of  Carstens,  Nellson  and  Altlnbilek  (1969)  where  T  -  3.56s.  At  this 
value  of  a  and  with  D  ■  0.585  mm,  curves  in  their  Figures  11  and  23  give 
A  *  27.0  cm  and  f^  *  .111  (with  their  f  replaced  by  3irf^). 


Figure  6.  Stress  coefficient  f  as  function  of  phase  0 
and  ripple  length  A,  for  natural  profiles. 
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IV.  DISCUSSION 

1 .  Conclusions  Regarding  f(e)  and  . 

Defined  by  the  magnitude,  shape  and  phase  of  their  peaks  and  valleys, 
the  character  of  the  curves  in  Figures  5,6,7  is  found  to  progress  contin¬ 
uously  from  one  to  the  next  as,  within  each  family,  U,  X  and  x/a  take  on 
successive  values.  This  continuity  supports  the  reliability  of  the  curves 
at  hand,  some  of  which  will  be  repeated  for  further  confirmation. 

Although  restricted  in  the  ranges  of  their  parameters,  the  curves  obtained 
permit  some  preliminary  conclusions  to  be  drawn.  It  is  to  be  remembered 
that  u  (0)  =  U  sin  0,  as  in  equation  (1). 

The  shapes  of  the  functions  f ( © )  in  Figures  5,6,7  differ  markedly  for 
different  values  of  the  parameters.  Therefore,  there  can  exist  no 
"universal"  function  f(e)  to  which  all  curves  can  be  reduced  by  a  selec¬ 
tion  of  scale.  In  particular,  in  no  case  has  f(0)  been  found  to  be  pro¬ 
portional  to  sin20.  Thus,  the  instantaneous  stress  is  not,  in  general, 
proportional  to  the  square  of  the  instantaneous  velocity. 

The  values  of  f^  listed  in  Table  2  are  found  to  vary  smoothly  with 
their  parameters  and,  so  far,  to  be  fairly  constant.  In  contrast,  f  , 
which  is  a^so  listed,  varies  more  abruptly  and  does  not  remain  propor¬ 
tional  to  f  .  Thus,  f  does  not  provide  a  reliable  measure  of  the  rate 
of  energy  dissipation. 

2.  No  Effect  of  Strain 


It  might  be  supposed  that  as  a  developing  profile  approaches  equili¬ 
brium,  or  as  a  strained  profile  approaches  an  unstrained  "natural"  con¬ 
dition,  the  stress  would  take  on  some  special  stationary  character, 
expressing,  in  some  form,  the  principle  of  lea_st  action.  However,  the 
curves  of  f(0)  in  Figure  7  and  the  values  of  f-  in  Table  2  vary  contin¬ 
uously  with  X/a,  with  no  discernible  special  character  or  extremum  at 
X/a  =  1.33  where  the  profile  is  unstrained  or  natural.  The  ripple  pro¬ 
files  of  Figure  7,  though  strained,  are  stable,  and  these  experiments 
have  not  examined  the  stress  on  unstable  profiles  during  times  of  growth 
or  change. 

3.  A  Simple  Model 

Certain  features  common  to  all  the  curves  of  f(0)  in  Figures  5,6,7 
suggest  a  general,  though  limited,  interpretation  in  terms  of  a  simplified 
physical  model.  These  features  are  three:  (1)  a  positive  stress  at 
0*0,  (2)  one  or  two  peaks  in  f(0)  in  the  range  0  4  9  <  tt/2,  and  (3)  a 
single  peak  in  f(0)  in  the  range  tt/2  <  0  <  tt.  It  is  first  recalled  that 
the  stress  on  a  rigid  flat  bed  precedes  0^(0)  by  a  shift  in  phase  which  is 
a  maximum,  tt/4,  for  a  smooth  bed  with  laminar  flow  and  stress  and 
decreases  as  the  bed  becomes  more  rough  and  the  stress  more  turbulent. 
Thus,  with  u  (0)  ■  U  sin  0  such  a  stress  would  have  a  maximum  in  the  range 
tt/A  <  0  <  n/?.  The  model  simply  combines  such  a  "flat-bed  stress"  on  the 
upstream  face  of  the  ripple  crest  with  the  stresses  of  vortex  action, 
which  is  essential  to  ripple  equilibrium. 
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The  process  of  vortex  formation  and  action  adopted  by  the  model  is 
familiar.  Following  flow  reversal,  tangential  stress  on  the  upstream 
ripple  face  generates  a  boundary-vorticity  layer  and  sets  sand  into 
motion,  skewing  the  crest  to  end  abruptly  in  a  steep  lee  slope  and 
sliding  and  spilling  into  the  trough.  The  sharp  curvature  of  the  crest 
enhances  separation  of  flow,  and,  as  it  is  carried  past,  the  boundary 
layer,  with  its  vorticity,  detaches  and  rolls  into  a  vortex,  still 
spilling  sand  into  the  trough.  The  center  of  the  vortex  is  slowed  and 
then  held  over  the  trough  by  reduced  pressure  on  the  lee  slope  of  the 
crest,  creating  a  fora  drag.  Around  the  time  of  the  succeeding  flow 
reversal,  the  spinning  vortex  sweeps  sand  out  of  the  trough  and  back  up 
or  over  the  crest,  thus  maintaining  the  profile. 

Returning  to  6  ■  0,  the  spin  of  the  vortex  formed  during  the  preced¬ 
ing  half  cycle  and  centered  over  the  trough  contributes  a  positive  stress, 
while  pressures  needed  to  hold  the  vortex  in  position  vanish  at  about  this 
time.  Thus,  f(0)  >  0.  Following  flow  reversal,  the  vortex,  decaying  but 
still  spinning,  is  carried  up  the  ripple  face  where  its  spin  hastens  and 
increases  the  stress.  Spin  decay  and  the  rapid  drop  off  in  spin  stress 
as  the  vortex  passes  the  ripple  crest  can  account  for  an  early  maximum  of 
the  spin  stress,  slightly  before  or  after  8  *  0,  and  preceding  the  maximum 
of  the  flat-bed  stress.  Thus,  depending  upon  the  relative  strengths  of 
these  two  maxima  and  their  separation  in  time,  f ( 0)  can  display  one  or  two 
peaks  in  the  range  0  <  8  <  it/2.  The  deceleration  of  the  forming  vortex 
(or,  more  accurately,  of  several  small  vortices  which  later  amalgamate) 
by  reduced  pressure  on  the  lee  slope  begins  with  flow  separation  soon  after 
flow  reversal,  but  the  associated  form  drag  reaches  a  single  peak  as  the 
vortex  approaches  full  development  sometime  in  the  range  tt/2  <  0  <  tt  ,  when 
the  flat-bed  stress  is  in  decline.  A  negative  back-flow  stress  on  the  lee 
ripple  slope  and  under  the  vortex  accompanies  this  fora  drag,  and  develops 
to  make  f (tt)  <  0  as  the  form  drag  vanishes. 

The  value  of  f_  from  Carstens,  Neilson  and  Altinbilek  (1969),  noted 
in  section  III,  4  above,  fits  well  with  values  listed  in  Table  2. 

However,  in  Figures  5,6,7  the  peaks  in  f(6)  preceding  and  nearest 
8  -  tt/2,  which  the  model  associates  with  the  flat-bed  stress,  are  an 
order  of  magnitude  larger  than  values  of  f  obtained  by  Kamphuis  (1975) 
at  the  same  values  of  a/D  and  aU/v.  At  least  to  some  extent,  this  large 
increase  may  be  a  result  of  the  movement  of  sand  grains  on  the  ripple 
face  (cf.  Lofquist  in  discussion  of  Vitale,  1979).  Similarly,  the  peaks 
in  f (0)  following  6  *  tt /2 ,  which  the  model  associates  with  fora  drag  and 
vortex  formation,  suggest  a  fora  drag  contribution  to  f  of,  roughly, 
one  third.  This  is  several  times  the  seven  percent  which  Tunstall  and 
Inman  (1975)  have  estimated  for  the  ratio  of  vortex  energy  to  the  total 
energy  dissipated  in  a  half  cycle.  Their  estimate  was  based  upon  a 
study  of  vortex  formation  behind  rounded  rigid  (lithlfied)  ripple  crests. 
Possibly,  the  relative  sharpness  of  a  naturally  skewing  crest  substan¬ 
tially  enhances  flow  separation  and  vortex  formation,  or,  the  vortex 
energy  is  only  a  part  of  the  fora  drag  contribution  to  f^. 
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4.  Limitations  and  a  Correction 


Whatever  its  merits,  the  model  is  severely  limited.  Processes 
which  interact  cannot  simply  be  added  as  flat-bed  and  vortex  stresses. 
Concentration  of  all  the  vorticity  generated  on  the  upstream  face  of 
a  crest  into  a  single  vortex  formed  in  its  lee  is  not  a  universal 
mechanism.  In  particular,  at  higher  rates  of  flow,  and  scour,  the 
ripple  crests  become  lower  and  closer  and  unable  to  form  and  hold 
vortices  large  enough  to  contain  all  the  detached  boundary  layer,  and 
flow  patterns  modify.  Such  appeared  to  be  the  case  with  exp.  28. 

When  large  quantities  of  sand  are  set  in  motion,  f(0)  as  calculated 
by  equation  (9)  itselt  Ic^js  precision  and  requires  the  addition  of  a 
correction,  Af(B),  .f,  at  any  0,  the  sand  in  motion  lowers  the  mean 
surface  by  Ah  (a  positive  quantity)  and  has  a  momentum  m  per  unit  area 
of  the  bed,  and  vfth  Q  *  bHU  sin  0,  a  rederivation  of  f?0)  produces  the 
correction  term 

Af (0)  =  —  cos  0  -  m  (15) 

a  pUz  s 

In  experiment  28,  Tigure  5  shows  f(8),  surprisingly,  to  be  negative 
in  the  neighborhood  of  0  =  tt/3.  Here,  Ah  cos  0  is  positive  while  m  is 
presumably  negL'L;ve,  as  sand  scoured  from  the  bed  during  the  preceding 
peak  in  f(0)  falls  back  to  impact  the  bed  and  lose  its  momentum.  Thus, 
Af(0)  is  presumably  positive  and,  while  its  magnitude  can  only  be  guessed 
at,  it  could  suffice  to  keep  the  corrected  f(0)  positive. 

5.  Future  Work 


The  experiments  described  here  are  part  of  a  larger  program  inter¬ 
rupted  by  the  timing  of  the  conference.  As  noted  above,  future  experi¬ 
ments  will  investigate  the  repeatability  of  present  results,  and  will 
extend  the  ranges  and  the  number  of  combinations  of  the  parameters.  One 
or  two  finer  sands  will  be  included.  Future  more  complete  analyses  will 
more  consistently  combine  parameters  into  dimensionless  ratios.  A 
completed  study  is  scheduled  for  publication  in  1981  by  the  Coastal 
Engineering  Research  Center,  Ft.  Belvoir,  Virginia. 
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